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1. INTRODUCTION 

Graphene is a two dimensional crystal composed of car¬ 
bon atoms which are packed in a honeycomb (hexagonal) 
lattice [1, 2]. It attracts considerable interest because of 
its unique electronic properties; most of them are related 
to existence of two conical points in the electron energy 
spectrum (Fermi points) and “massless Dirac fermion” 
character of charge carriers with energy and momentum 
close to the Fermi points [3-7]. It results in numerous 
quantum relativistic phenomena such as Klein tunnel¬ 
ing, minimal conductivity through evanescent waves, rel¬ 
ativistic collapse at a supercritical charge, etc., estab¬ 
lishing an interesting and fruitful relation between funda¬ 
mental physics and materials science [8-13]. The effective 
“velocity of light” (Fermi velocity) for the Dirac fermions 
in graphene is relatively small, vf ~ c/300, and the in¬ 
teraction between the quasiparticles in graphene can be 
approximated by the instantaneous Coulomb potential 
with the effective coupling constant * 1 a e ff ~ a(c/vF ) ~ 
300/137 ~ 2. This interaction is therefore quite strong, 
which results in a rich variety of phenomena [12]. Within 
the Dirac model, before the experimental discovery of 
graphene, it was shown that the long-range character of 
Coulomb interaction results in a renormalization of the 
Fermi velocity which is divergent at zero temperature and 
zero doping leading to a non-Fermi-liquid behaviour [14]; 
this prediction has been recently confirmed experimen¬ 
tally [ 5, 16]. 

At the same time, Dirac model gives us many-body 
renormalization of electronic properties only for small 
coupling constant and only with the logarithmic accu¬ 
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racy. The higher-order terms were considered in Refs. 
[17, 18] but still within the Dirac model. To calcu¬ 
late quantitatively correctly these properties one needs 
to work with the lattice model and with a realistic po¬ 
tential of electron-electron interaction V xy taking into ac¬ 
count its screening by the er-electrons. The correspond¬ 
ing first-principle results [19] can be parametrized by the 
phenomenological potential. This modification of the in¬ 
teraction potential at small scales in comparison with 
the bare Coulomb potential was proved to affect some 
graphene properties significantly (for instance, the phase 
diagram of graphene [20]). 

The authors of papers [20-22] carried out the Monte- 
Carlo study of graphene properties based on the tight- 
binding Hamiltonian without the expansion near the 
Fermi points. Within this approach one can introduce 
an arbitrary phenomenological potential V xy . Using the 
tight-binding Hamiltonian on the hexagonal lattice, we 
are going to build the perturbation theory in V xy . We 
believe that the theory built in this way is important 
and interesting for the following reasons. 

The theory based on the tight-binding Hamiltonian has 
more common features with real graphene physics than 
the effective theory based on the expansion in the vicinity 
of Dirac points. For instance, the tight-binding Hamilto¬ 
nian “remembers” graphene properties such as geometry 
of hexagonal lattice or the natural energy scale such as n- 
bandwidth, which are absent in the effective Dirac theory. 
Moreover, one can include a phenomenological potential 
in this theory, which is closer to the real graphene physics 
than the bare Coulomb potential. 

In addition, the theory based on the tight-binding 
Hamiltonian and the phenomenological interaction V xy 
can be easily improved. For instance, one can study 
the effects appearing due to the inclusion of the next-to- 
the-nearest-neighbour hopping or nonzero chemical po¬ 
tential. Note that such study cannot be carried out 
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within the Monte-Carlo simulation because of the well 
known sign problem. Finally, if the electron properties 
of other nanomaterials are formulated in terms of the 
tight-binding Hamiltonian and the phenomenological po¬ 
tential, it makes no difficulties to apply the results of this 
paper to study these materials. 

Lattice simulation of graphene were proved to be a very 
efficient and quickly developing approach for studying the 
properties of graphene [20-25]. An important feature of 
all these simulations is that they are conducted at the fi¬ 
nite lattice, at the finite temperature and with the finite 
discretization errors. In order to check the lattice results 
and estimate the discretization uncertainty in the weak 
coupling region it is very useful to develop the perturba¬ 
tion theory which accounts all these effects. 

Strictly speaking, the theory with the arbitrary phe¬ 
nomenological potential is not renormalizable since it 
contains four-fermion terms. However, lattice formula¬ 
tion provides the ultraviolet (spacing between carbon 
atoms) and the infrared (the finite size of the lattice) 
regulators. For this reason the theory on the hexagonal 
lattice is well defined. 

Finally, one should mention that the interaction in 
graphene is strong, so the application of the perturbation 
theory is questionable. However, we believe that even at 
the one-loop level one can study some important physical 
effects. One can also expect that the perturbation theory 
built in this paper can well describe graphene many-body 
effects similarly to the RPA based on effective theory of 
graphene[18]. 

This paper is organized as follows. In the next section 
we built the perturbation theory which is based on the 
tight binding Hamiltonian and arbitrary interaction po¬ 
tential between electrons. In the section 3 we calculate 
one-loop corrections to the electron propagator. The sec¬ 
tion 4 is devoted to the calculation of one-loop corrections 
to the interaction potential. Finally, in the last section 
we discuss and summarize the results of this paper. 
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FIG. 1: L x x L y = 6 x 6 graphene sheet. The sites 
belonging to the sublattices A and B are shown as 
rectangles and circles respectively. 


£i 0... L x 1, £2 0... Ly 1, so 


— \/3a£i + 
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(1) 

— ^ aS s ,B■ 

(2) 


The torus topology implies the following identification of 
£ coordinates: 


(Ci + L x , £2) —t (£ 1 , £2)) (£i, £2 + L y ) — > (£i + Lyj 2 , £2). 

(3) 

Every site of the A sublattice is connected to three 
sites of the sublattice B. The vectors pb in £-coordinates, 
connect (a,£) to its neighbours (/3,£ + pi): 

Pi = (0,0), P 2 = (-1,1), P3 = (- 1 , 0 ). (4) 


In x-coordinates these pi vectors read: 


2. PARTITION FUNCTION AND FEYNMAN 
RULES 


Pi = (0,0), p 2 = ” 


h , p 3 = (-\/3a,o) . (5) 


2.1. Geometry 

We consider a hexagonal lattice with the torus topol¬ 
ogy. The example of such lattice, which consists of 
L x x L y = 6 x 6 hexagons, is shown in Fig.l. 

The hexagonal lattice is the composition of two tri¬ 
angular sublattices A and B. The sites belonging to 
the sublattices A and B are shown as rectangles and 
circles respectively. The Cartesian coordinates ( x,y ) of 
any lattice site can be parametrised by three numbers 
(s, £i, £ 2 ), where s = A/B is the sublattice index and 


2.2. Tight-binding Hamiltonian with interactions 

The electronic properties of graphene can be described 
by the tight-binding Hamiltonian: 

^tb = ^ ' y 1 (^t,x^,y + Ca : yCtT, X ) , (6) 

<y=t,i {xy) 

where the summation is kept over the neighbouring 
graphene lattice sites x and y, t ~ 2.7 eV is the hop¬ 
ping parameter. Operators c' a x , c„ tX create and annihi¬ 
late electron with spin a =f, 4- at the lattice site x. Note 









































3 


that one can also include the next-to-nearest neighbours 
hopping to the Hamiltonian (6), but for simplicity we 
restrict our consideration to the nearest neighbours. 

We choose the vacuum state that satisfies the following 
conditions: 


Cf,x |0> = 0, c} iX |0> = 0, (7) 

so there is an electron with spin cr =1 at every lattice 
site and no electrons with spin er ='|\ As our main goal 
is to calculate the partition function Z , which contains 
summation over all states, the specific choice of |0) will 
not affect any physical results. 

It is convenient to rewrite the Hamiltonian in terms of 
creation and annihilation operators for ’’particles” and 
’’holes”, which read: 

b x — C_|-,x, bx — Ac_ x . (8) 

The plus sign is taken for x € A and the minus sign for 
x € -B, where A and B are the triangular sublattices la¬ 
bels. After operators redefinition the ground state satis¬ 
fies a x |0) = b x |0) = 0. Thereby, we interpret the absence 
of a valence electron as a positively charged ” hole” and an 
additional electron as a negatively charged ” particle”. In 
terms of these operators the tight-binding Hamiltonian 
is: 


H =—t (a) x a y + P y b x + H.c.^j . (9) 

<xy> 

The charge operator now reads: 

Qx ® X Q'X ®x^x- ( 10 ) 

It easy to check that q x |0) = 0, which means that this 
vacuum is electrically neutral. 

In some applications [20, 21] the mass term is intro¬ 
duced: 

H m = '^{±mc 2 )a\a x + '^2{±mc 2 )b\b x , (11) 

X X 

here the plus sign is taken for sublattice A and the minus 
sign — for sablattice B . This term explicitly breaks the 
symmetry between two sublattices. 

In order to study the action of the chemical potential 
on graphene properties we introduce the term 

Hy, = M {a\a x - b\. b x ). (12) 

X 

It is known that the interaction between electrons plays 
an important role. This part of the Hamiltonian has the 
form: 

-^int = 9 VbvQbQv (13) 

x,y 

The Coulomb potential V is often used to describe the 
interaction. However, it was shown in [19], that the real 


potential V xy dramatically deviates from the Coulomb 
law at small distances, which affects physical properties 

[ 20 ]. 

The aim of this paper is to study the properties of the 
electronic system, described by the Hamiltonian: 

H = —t ^2 \®i®v + x + H.c.j + y^(±mc 2 )al.a x + 

<xy> x 

A 'y ^ (zfcrnc )b x 6 x + /i y ] (ci\.a x — b^bx) + ^ ^ ) VxyQxQyi 

x x x,y 

(14) 

where V xy is an arbitrary phenomenological interaction, 
treated as the perturbation. 


2.3. Electronic spectrum of graphene without 
interaction 

If the interaction is neglected 2 , one can easily find 
Hamiltonian spectrum and its eigenfunctions, which can 
be written as: 

(15) 

where £ = ±1 is an additional label representing par¬ 
ticles (C = +1) and antiparticles (£ = —1). The first 
vector component corresponds to the A sublattice, simi¬ 
larly the second one corresponds to the B sublattice. The 
Brillouin zone momentums k are: 

k,=^,k,= 2 ; m » (i6) 

v'Sai, " ML,/2) ' ’ 

Indices m x = 0,..., L x — 1, m y = 0,..., L y — 1 give the 
full set of the eigenfunctions for the torus topology of the 
graphene sheet. 

The energy spectrum of the Hamiltonian is then: 

E^(k) = C,\J m 2 c A + \(fi(k)\ 2 , where ip(k) = t e*^. 

6=1 

(17) 


The wave function components are 

4(fc) = 


I E(k ) + £ me 2 


2 E(k)L x L y 


<4(£) = - 


Ce —i mgy (S) E(k)-(mc 2 ' 
V 2 E{k)L x L y 


There exist two special momenta kp, such that 
y(k F ) = 0 and E(kp ) = ±rac 2 , called the Dirac points. 


2 In this subsection we consider graphene with zero chemical po¬ 
tential 
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There are two different Dirac positioned at: (m x ,m y ) = 
(2L x /3,0) and (L x /3,L y /2). 

For the massless fermions (m = 0), in the vicinity 
of the Dirac points the energy spectrum of the fermion 
quasiparticles is: 

E(k) = Vp\k — kp |, 

v °f = ^ = 0.87-10 6 (19) 

2 344 s 

The linear spectrum of fermion excitations plays the cen¬ 
tral role in the effective theory of graphene. 

2.4. Partition function 

Let us consider a graphene sheet L x x L y at the tem¬ 
perature T. To write down the partition function for this 
graphene sample, one should handle the 4-fermion inter¬ 
action operator which is contained in the Hamiltonian 
(14). So, before going to the partition function the in¬ 
teraction term of the Hamiltonian should be decomposed 
using the Hubbard-Stratonovich transformation: 



-►-electron 

- M -hole 

’’photon” 

FIG. 2: Graphical representation of the propagators of 
different particles. 


In the last expression we introduced the Hubbard field 
p, which carries the interaction. Note also that we have 
omitted the determinant of the V which is not important 
for our calculations. 

Now it causes no difficulties to write down the partition 
function in terms of the path integral. To this end, we 
divide the Euclidean time Tp £ [0,/3], /? = 1/T into N T 
parts with the size St = 8/N T . We introduce the electron 
fields (fj, rj), the hole fields (ip, ip) and the Hubbard field p 
(which is called “photon” below) at each site of the lattice 
L x x L y x f3. The boundary conditions in Euclidean time 
direction are periodic for “photon” p and antiperiodic for 
the fermion fields (fj, rf), (ip, ip). The boundary conditions 
for the spatial directions are periodic for all the fields. 

Finally, the partition function for the model can be 
written as [20, 21]: 


Z = j 'Dp'DvfDrj'Dip'Dipe a ' x ' v a ' u: ' * * v , 


( 21 ) 


where the matrix M,i ll2a (x,y) reads (here l\ = A,B and Z 2 = A,B are sublattice indices): 


Mi lt i 2 (x,y) = ( 5 St y6i lt i 2 6 x o ty o - e l5ripi ^ x) e Sr ^ 5s^5i lM 5 x o +ST , y o +tSr e' 5r/i 


e iSTip A (x)S h B ^2 bxO+5T,y°bx+pi,y+ 


+e lSTIpB ^5l ltB 8l 2 ,A ^2 $E°+<S T,y°dx-Pi,y 


i=l 


+ me 2 St e ST » [e lSrVA{x) 5 hA 5 l2A 5 s ^5 x0+ST , y o - e iST ^6 htB 6 l2tB S s , y 5 x0+ ^ y o 

( 22 ) 


St 


Sem(<p) = — ^2p(x) V xy 1 p(y). 


(23) 


To get the expression for the fermion operator M one 
should take the formula for the M and carry out the 
substitution p -4 — p , p —> —p. Note that in the for¬ 
mula (23) the summation is taken over all the coordinates 
(£1 ,£ 2 ) and over the sublattice indices, so the interaction 


potential V xy is assumed to be a matrix 

_ (V AA (x,y) V AB (x,y)\ 

v *y-\yBA {x y) V BB {x ^ y) j ■ 


(24) 


In the limit of St —> 0 expression (21) corresponds to 
the partition function of graphene sheet L x x L y at the 
temperature T. In some applications, for instance in the 
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FIG. 3: The interaction vertex of electron with one 
“photon” field 

Monte-Carlo simulation of graphene [20, 21], one uses 
small but finite step St. Some of the formulas below are 
written for the finite St to have the possibility to estimate 
the discretization uncertainty. 


FIG. 4: The interaction vertex of electron with two “pho¬ 
ton” fields 

2. The hole propagator ( rj(y)fj(x )) can be obtained 
from the electron propagator with the substitution 

3. The propagator of the “photon” field <p 


2.5. Propagators 

Using the expression for the partition function (21), 
one can write down free propagators for the correspond¬ 
ing fields (see Fig. 2). 


(ip{x,x°)ip(y,y 0 )) = ^-S x o iy oV(x,y) = 

1 y r ik (x -y) (V AA (k) V AB (k)\ 
St L T L x L y \V BA (k ) V BB {k)J ' 


1. The electron propagator 


(v(y)v(x)) = M^(x,y) = 

1 ^ e ik(x-y) 

L t L x Ly (1 _ e ik° St & 6t [i\‘2 _ § T 2 £2 ( fc\ e 2ik° Sr £ 25t n * 
k° ,k 

(l — (1 + mc 2 (5r)e* fe ST e 2Sr t 1 ip(k)STe lk ° ST 

V ip*(k)STe ik ° ST 1-(1 -mc 2 ST)e ik ° 5T e STfl 

(25) 


The potential is instantaneous and acts in one Eu¬ 
clidean time slice. 

Using the expressions for the partition function, for the 
electron and hole propagators one can show that in the 
limit St —>• 0 the charge of graphene sheet (Q) is 


(Q) = - E ( Tr ( M -) - Tr fc 1 )) = 

x ' ' 

= 2 Vf-^- J- - 

^ y e p(E(k)-») + l e f>(E(k)+n) + i 


(26) 


where i p{k) was introduced in (17). 

Note that the matrix M^ 1 {x, y) is the 2x2 matrix 
in the ’’sublattice space”, 


M 0 1 (.r. y) 


(M^ A {x,y) M^ B (x,y)\ 
\M BA {x, y) M BB {x, y)j 


its elements correspond to the propagation between 
one sublattice (AA and BB) and different sublat¬ 
tices (AB and BA). 


as it should be. 


2.6. Vertices 


The partition function (21) takes into account the in¬ 
teractions between electrons/holes and “photons”. The 
interaction vertices of “photons” and electrons can be de¬ 
rived with an arbitrary accuracy through the expansion 
of the expression for the fermion operator: 


SMi ltl2 (x,y) =(e i5Tt *iM - l) e* T " -S^S llM S x o + ^ y o + tSr 


^l 2 ,B ^ ^ ^x 0 -\-Sr,y°^x-\-pi,y 
i =1 


^ ^ fix°-\-8T,y°&x—pi,y 


-|- TYIC St ,aSi 2 ,ASx,ySx°-\-Sr,y° &li ,bSi 2 ,bSx, ySx°-\-8T,y°\ 


(27) 
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► 



FIG. 5: Diagram that contributes to the self energy 

e i(p) 


FIG. 6: Fock diagram that contributes to the self 
energy E 2 (p) 


in powers of St. The simplest interaction vertex is At the leading order approximation there are two dia- 
shown in Fig. 3. Using the formula (27), one can show grams shown in Fig. 5, 6 that contribute to E(p): 
that 


at the leading order approximation in <5r-expansion 
this vertex can be written as follows 


s i (p) = +2 6t 


1 , fV Q0 0 


00 

0 U 00 


P ip° St 


V ( 3 ) =(-iSt ) Y (S s ^Si u i 2 S x o + s T ^yo)ip h {x)ip h {x)ipi 2 {y) S 2 (p) = -Si (p) + £ f T V tanh ( I x 

X V 

(28) 


2T 


where l\. Z 2 are the sublattice indices. 

The next vertex, which is required for the subsequent 
analysis and describes the interaction between two elec¬ 
trons and two “photons”, is shown in Fig. 4. 

At the leading order approximation it can be written 
as: 

<57~2 

y (A) = ~tt Y ( S S,ySl lt l 2 S x o +STl yo)'lf)l 1 (x)ip , f 1 (x)lpl 2 (y) 


h,h,x,y 


(29) 




The inverse of the free electron propagator Mo(p) can 
be written in the following form (see formula (25)): 


M 0 (p) = 


\ _ e ip°Sr 

0 


0 

l _ g ip°Sr 


+ e ip ° ST Sr 


me 2 (p(p) 
ip*(p) —me 2 

(31) 


The final expression for E(p) is 


There are also additional vertices, coupling electrons 
to 3,4..., n,... “photons” (because the ip field stands 
in the exponent). However, they are suppressed by the 
additional factors St and give no contribution to the final 
answer in the limit St -a 0. Note also that in this section 
we presented only vertices with electrons. The vertices 
with holes can be found similarly. 


£(„) = E,(p) + S,(P) = PPP J>„h l E & 


2 L X Ly 


2 T 


mv^[p-ky 


V* (fc) \fBA 
E(k) 


(p~ k ) ^y BB {p-k) l 


(32) 


3. ONE-LOOP CORRECTIONS TO THE 
ELECTRON PROPAGATOR 


First we are going to consider one-loop corrections at 
zero chemical potential and nonzero temperature. One- 
loop corrections to the electron propagator can be ex¬ 
pressed in terms of the self-energy function E(p): 


M~\p) 


1 

M 0 (p) - £(p)' 


(30) 


From the formula (32) one can see that one-loop cor¬ 
rections are reduced to the renormalization of the mass 
m and the function (fi(p). The renormalized mass now 
depends on the point in the Brillouin zone and has the 
form 


m R (p) = m + 


2 L x Ly 


Y^ tanh 


m 

2 T 


V AA (p- k). 


E(k ) 


(33) 
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(b) The spectrum profile is drawn on 

(a) Free and renormalized energy E(k) profile. the polygon ABCD. 

FIG. 7: Energy spectrum profile of electrons with Coulomb, screened Coulomb [19] interactions and without 
interactions. The calculation was carried out at T = 0.1 eV at infinitely large lattice. 


The expression for the renormalized function <p R (p) is 
‘ P R (P) = <P(p) d 

- k v 

-V AB (p-k). (34) 


PW -TrABf 


1 ^tanhf ® 1 x 


2 T 


E{k) 


Thus one-loop corrections conserve the form of the prop¬ 
agator (25) and substitute the free mass m and the func¬ 
tion ip(p) with the renormalized expressions (33), (34). 
From this one can conclude that the energy spectrum 
of the quasiparticles at the one-loop approximation is 
E 2 = ( m R c 2 ) 2 + (ip R ) 2 ■ I n order to estimate the size of 
one-loop corrections in Fig. 7 we plot the energy spec¬ 
trum profile of electron with Coulomb, screened Coulomb 
[19] interactions and without interactions. The calcula¬ 
tion was carried out at T = 0.1 eV at infinitely large 
lattice. 

The formulas (33) and (34) can be used to reproduce 
well-known results of the effective theory of graphene and 
generalize them to the case of nonzero temperature. To 
this end, we consider large lattice L x , L y —> oo with the 
Coulomb interactions between electrons near the Fermi 
point. Then the formulas (33) and (34) can be written 
as 


m R (p = p F ) = 


, , 1 a 

= m 1 + - - 


2 (v F /c) 


, , 1 a 

= v F 1 + -- 


4 {y F /c) 


log 


log 


v F A 
2 cT 

v R - 
V F — 

v F A 

2 cT 


+ 7 - log 7 r /4 + 0(A 1 ) 

+ 7 - log 7 t /4 + 0(A~ 1 ) 
(35) ' 


where A is the ultraviolet cut-off, 7 ~ 0.577... is the Eu¬ 
ler’s constant. These formulas are in agreement with the 
predictions of effective theory of graphene at the leading 
logarithmic accuracy [ 4, 26]. 

During the derivation of the formulas (35) we carried 
out the integration over the two-dimensional sphere of 
radius A. It was also assumed that me 2 CT, A is much 
larger than all other energy scales involved. Note that the 
theory is regularized by the temperature in the infrared 
region. Note also that the ratio of the ultraviolet and the 
infrared cut-offs A/T under the logarithm is multiplied 
by the Vf/c that considerably reduces the total renormal¬ 
ization of the Fermi velocity. 

In order to study the finite density effects we consider 
one-loop corrections at zero temperature and nonzero 
chemical potential. In addition to the diagrams shown 
in Figures 5, 6 there is a contribution of the interaction 
of the electron propagator with vacuum represented by 
the diagram shown in Fig. 8 . As before, one-loop cor¬ 
rections at nonzero p lead to the renormalization of the 
parameters of free propagator preserving its structure. 
In particular, the mass and the function ip at one-loop 
approximation can be written as follows 


(p) = m+ m TAt\V AA {p- k) 6(E(k) - p), 


IP ($>) = P(P) + 


2L. r L 

1 


x^y E(k) 


j2^-v AB (p-k) e{E{k)-p). 


2 L x L y E(k) 

k 


(36) 


In the limit of large lattice L x , L y —> 00 , Coulomb inter¬ 
action and the linear electron spectrum near the Fermi 





















FIG. 8: Hartree diagram that 
contributes to the one-loop 
electron propagator at p ^ 0. 


FIG. 9: Empty-loop diagram that 
contributes to the renormalization 
of the interaction potential. 


point, the formulas (36) can be written as 


m H (p = pp) = m 1 + 


1 a 


2 (v F /c) 


Vp = Vp 1 + - 


1 a 


4( v F /c ) 


log 


log 

v f A 

cp 


v f A 

cp 


(37) 


Similarly to the case of nonzero temperature chemical po¬ 
tential p plays the role of the infrared cut-off. It should 
be noted now that there are four scales that can play 
the role of the infrared cut-off: the fermion mass to, the 
chemical potential, the temperature and the inverse lat¬ 
tice size. We believe that in the infrared limit the theory 
is regularized by the largest of these scales. Note that 
the formulas (35), (36) take these effects into account 
exactly. 

Note also that the formulas (35) can be written in the 
form (37) if instead of A one uses new ultraviolet cut off 
A' = 2e 7 /7r - A ~ 1.13 • A and substitute the chemical 
potential by the temperature. 

For this reason at temperature equal to chemical 
potential the renormalization due to temperature effects 
is a little bit larger than the renormalization due to 
nonzero chemical potential. In this consideration we did 
not take the Debye screening of the interaction potential 
at large distances and the screening of the Coulomb 
potential at small distances into account, that will be 


done in the last section. 


4. ONE-LOOP CORRECTIONS TO “PHOTON” 
PROPAGATOR AND DIELECTRIC 
PERMITTIVITY OF GRAPHENE 

In this section we calculate one-loop correction to the 
instantaneous “photon” propagator: ((p(x,x°)ip(y,y 0 )'). 

One-loop corrections to the propagator in the momen¬ 
tum space can be expressed in terms of the polarisation 

, ^ (V AA V AB \ 

operator V = I I , as 

V R (k) = -kr x V(k) 

m 

e(k) = l-V(k)xV(k), (38) 

where V R (k),V(k) are renormalized and tree level 
Fourier transforms of the potential, which are the 2x2 
matrices in the sublattice space. In the limit of dr —> 0 
and one-loop approximation the only diagram shown in 
Fig. 9 contributes. 

The expression for the polarization operator at one- 
loop approximation can be written in the following form: 


V AA {k)=V MB {k) = 


L X Ly 


E 


1 

E 2 (p- k) - E 2 (p) 


m 2 c 4 + E 2 (p) E(p) m 2 c 4 + E 2 (p—k) E(p — k ) 

tanh- = -tanh- 

27 1 E(p— k ) 2T 


m 


V™d) = v°\k) = jT- Y. 


m 


i 


E 2 (p-k) - E 2 (p)\E(p) 2 T E{p-k) 


tanh 


E(p— k ) 
2T 


(39) 


Using formulas (38), (39) one can show that at large dis¬ 
tances and at small temperature the expression for the 
interaction potentials for all the sublattice indices take 
the form 


V R (k ) 


2irac 1 
£ o (|fc| + todc) ’ 


(40) 


where 


eo = 1 + 


2 (vp/cY 
8a log2. 


m.D = 2 

eo v F 


-T. 


(41) 


Now one may see that at zero temperature the interac- 
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tion potential is the Coulomb potential screened by the 
dielectric permittivity eo, which is in agreement with the 
RPA result [13]. At nonzero temperature there is a two- 
dimensional Debye screening of the Coulomb potential 
with the Debye mass mo, which agrees with the results 
of the paper [27]. 

Similarly one can study the question how the nonzero 
density acts on the dielectric permittivity of graphene. 
To get the expression for the polarization operator in 
this case one can use the formulas (39) and carry out the 
following substitution tanh(.E/2T) —> 9(E — /i). It is not 
difficult to find out that at large distance the expression 
for the interaction potential has the form (40) with e 
given by the equation (41) and the Debye mass 

m D = -^\n\. (42) 

e 0 V F 

The last expression agrees with the RPA result [13]. 

At the end of this section we plot the dielectric permit¬ 
tivity of graphene e obtained from formulas (38), (39) as 
a function of distance r in units of the lattice spacing a 
(Fig. 11). The calculation was carried out for suspended 
graphene with the interaction potential from paper [20] 
and for different external conditions: T = 0, /i = 0, 
T = 0, n = 0.026 eV (n - 6.5 • 10 10 cm" 2 ) and T = 0.026 
eV, /x = 0. In order to compare our results with RPA 
we plot the dielectric permittivity given by formula (41). 
For small distance we compared our results with the re¬ 
sults of Monte-Carlo simulation [28] . It is seen that the 
Monte-Carlo results are in a good agreement with one- 
loop results. 

From Fig. 11 it is seen that at zero temperature the e 
starting from the value ~ 2 at r = 0 approaches the value 
eo (41) already for r/a > 5. At nonzero temperature 
and chemical potential the larger the r the larger the e. 
This behaviour can be attributed to the Debye screening 
effect. 


5. DISCUSSION AND CONCLUSION 

This paper is devoted to the perturbation theory which 
can be used for studying the properties of graphene at fi¬ 
nite temperature and nonzero chemical potential. This 
perturbation theory is based on the tight-binding Hamil¬ 
tonian on hexagonal lattice and arbitrary interaction po¬ 
tential between electrons, which is considered as a per¬ 
turbation. We built the partition function for this theory, 
derived Feynman rules and expressions for free propaga¬ 
tors. 

As an example of the application we calculated one- 
loop corrections to the electron propagator. It was shown 
that one-loop corrections lead to the renormalization 
of the bare mass and the function ip(k) conserving the 
structure of the propagator. Using this result we calcu¬ 
lated one-loop energy spectrum of electrons, renormal¬ 
ized fermion mass and Fermi velocity. 


In order to estimate the value of the renormalization, 
in Fig. 10 we plotted the renormalization factor for 
the Fermi velocity as a function of the temperature for 
graphene on hBN for the Coulomb and screened at small 
distances Coulomb interactions [19]. During the calcu¬ 
lation we took the screening of the potentials by the di¬ 
electric permittivity (38) into account. The results can 
be well described by the formula 

Vf(X)/v°f = (l + A-log 

where for the Coulomb interaction A = 0.096, A = 3.2 eV 
and for the screened Coulomb interaction A = 0.093, A = 
2.4 eV. The effective ultraviolet cut-off A is sensitive to 
the values of the potential at small distances, contrary 
to the coefficient A. Note also that that the coefficient 
A is well described by the expression A ss ac/(4v F e') = 
0.094, e '= (e BAr + l)/2 + e 0 - 1 = 6.7. 

In addition we considered one-loop corrections to the 
electron propagator at zero temperature and nonzero 
chemical potential. It was shown that similarly to the 
nonzero temperature case one-loop corrections conserve 
the structure of the propagator leading to the renormal¬ 
ization of already existing parameters. In particular, 
we obtained the formulas for the renormalization of the 
Fermi velocity and the fermion mass in this case. 

In paper [16] the renormalization of the Fermi velocity 
was studied through the measurement of quantum capac¬ 
itance at nonzero chemical potential. In this study the 
graphene layer was placed inside hBN, which reduced the 
strength of the interaction between electrons. This allows 
us to expect that perturbation theory works well in this 
case and we can compare our results with the the result 
of [16]. The results of the measurements of the Fermi 
velocity can be well fitted by the formula 



v F (n) = v f (/j, 0 ) 



1 a 
4 e(v F /c) 



= v F (no) 


^1 + A log 


Mo 

M 


(44) 


where Mf(mo) = 0.85 x 10 6 to/s, /jl o = 3.2eV (no = 
10 15 cm -2 ), e ~ 8 and A = 0.081. Note that the original 
formula for the Fermi velocity from [16] depends on the 
density n. In (44) we turn the dependence on n to the 
dependence on chemical potential. 

In the calculation we used the formulas (36) with 
the interaction potential screened by the dielectric per¬ 
mittivity (38). The interaction potential at small dis¬ 
tances was taken from the paper [19] and divided by 
€bn — 4.5. At large distances we took the Coulomb 
potential screened by cbn- The results of the calcula¬ 
tions can be well fitted by formula (44) with the pa¬ 
rameters: fj ,o = 2.9 eV, A = 0.072. This values are in 
reasonable agreement with paper [16]. In addition we 
carried out the calculation of the parameters fio,A for 
the Coulomb potential screened by cbn- Our result is 
fj ,o = 5.2 eV, A = 0.077. Again we see that the value of 
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FIG. 10: The renormalization factor for the Fermi 
velocity as a function of the temperature for graphene on 
hBN for Coulomb and screened at small distances 
Coulomb interactions [19]. 


10 
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6 

yj 
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2 


0 


FIG. 11: The dielectric permittivity of suspended 
graphene e as a function of distance r in units of distance 
between the neighbouring carbon atoms a for different 
external conditions. The dotted line corresponds to the 
RPA expression (41). The points shown as red triangles 
are the results of Monte-Carlo simulation. 


MC simulations i— *— 

effective theory. 

T = 0, /x = 0 • 

T = 0, fi = 26 meV x 
T = 26 meV, fi = 0 * 


:. 1.1 


if" 


-40“ 


20 


50 


/a 


the constant /uo is sensitive to the values of the potential 
at small distances. 

It is also interesting to compare the renormalization 
of the Fermi velocity due to nonzero /r and nonzero tem¬ 
perature. To this end we calculated nonzero temperature 
renormalization of the graphene layer placed inside the 
hBN for the Coulomb and screened at small distances 
Coulomb interactions. The results can be described by 
formula (43) with the parameters: for the Coulomb in¬ 
teraction A = 0.075, A = 4.5 eV and for the screened 
Coulomb interaction A = 0.073, A = 3.05 eV. 

The other example of application of the perturbation 
theory is the calculation of one-loop corrections to the 
interaction potential done in the previous section. We 
derived the one-loop expression for the dielectric permit¬ 
tivity at nonzero temperature, nonzero chemical poten¬ 
tial and the arbitrary interaction potential. 

It is well known that the electrons in graphene form a 
strongly interacting system. So it is reasonable to con¬ 
sider the question how our results are affected by the 
higher order corrections. The authors of paper [ .8] con¬ 
sidered the next-to-leading order (NLO) corrections to 
the Fermi velocity within effective theory of graphene. 
The main result of this paper in the statement that if 
one expands the Fermi velocity renormalization in the 
one-loop RPA potential instead of the usual Coulomb 


potential, the NLO corrections to the leading order(LO) 
result (37) turn out to be small. This allows us to expect 
a good accuracy of the formulas (33), (34), (36) with one- 
loop potential given by the formulas (38), (39) even for 
suspended graphene. 

The authors of the paper [17] considered the NLO cor¬ 
rections to the polarization operator within the effective 
theory of graphene. The NLO value of the dielectric per¬ 
mittivity for suspended graphene is approximately by 30 
% smaller than the LO result, which is not very large. 
Moreover, the Monte-Carlo results [2£] tell us that the 
higher order corrections to the LO result can be even 
smaller than 30 %. For this reason one can expect that 
the accuracy of formulas (38), (39) is rather good. 
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